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Abstract. We perform an exhaustive study of the role neutrino telescopes could 
0^ i play in the discovery and exploration of supersymmetric extensions of the Standard 

Model with a long-lived stau next-to-lightcst superparticle. These staus are produced 

in pairs by cosmic neutrino interactions in the Earth matter. We show that the 

^ I background of stau events to the standard muon signal is negligible and plays no 

OO . role in the determination of the cosmic neutrino flux. On the other hand, one can 

expect up to 50 pair events per year in a cubic kilometer detector such as IceCube, if 
•^^ ' the superpartner mass spectrum and the high-energy cosmic neutrino flux are close to 

^^ . experimental bounds. 

o 

rS^ ' E-mail: markus.ahlers@desy.de, joern.kersten@desy.de, 



andreas . ringwald@desy . de 



Long-Lived Staus at Neutrino Telescopes 2 

1. Introduction 

Supersymmetry (SUSY) is currently the most popular extension of the Standard Model 
(SM) and predicts a superpartner for each SM particle. In most models, i?-parity is 
conserved, so that the lightest superpartner (LSP) is stable. This makes it an excellent 
candidate for the dark matter. Most studies assume that this particle is a neutralino, 
which interacts weakly and may be observed therefore in direct dark matter searches. 
However, if SUSY is extended to include gravity, there is an alternative LSP candidate, 
the gravitino. As it is the superpartner of the graviton, it takes part only in the 
gravitational interaction. Therefore, the decay of the next-to-lightest SUSY particle 
(NLSP) is highly suppressed.! If it is a charged particle such as the stau, tr, the 
superpartner of the right-handed tau, it can possibly be collected in collider experiments. 
Observation of the stau decays could then lead to an indirect discovery of the gravitino. 
This exciting possibility has attracted considerable interest recently [SI 111 El El C] • 

Within this scenario, we elaborate in this paper on an alternative experimental 
approach in which the high energy necessary to produce SUSY particles is not provided 
by man-made accelerators but by nature in the form of high energy cosmic neutrinos. 
If the stau NLSPs resulting from neutrino-nucleon interactions inside the Earth are 
sufficiently long-lived, they can be detected in large ice or water Cherenkov neutrino 
telescopes, as pointed out in . These observatories are designed to measure the flux of 
high energy cosmic neutrinos via Cherenkov photons emitted by upward-going muons 
and cascades, which are produced by weak interactions in the Earth. High energy muons 
are visible if they originate up to a few tens of kilometers outside the detector. This 
defines the effective detection volume, which is limited by the energy loss of muons 
in matter. For a stau this effective detection volume increases dramatically due to 
the much smaller energy loss in matter [8^. This might compensate the suppression of 
the production cross section of SUSY particles compared to the one of the SM weak 
interactions. Moreover, interactions of cosmic neutrinos with nucleons inside the Earth 
will always produce pairs of staus, which appear as nearly parallel muon-like tracks in 
the detector due to their large boost-factor [Z| . This is expected to be in contrast to the 
SM processes, which lead to muon pairs only in rare cases. 

The goal of this paper is to perform an exhaustive study of the role which neutrino 
telescopes could play in the discovery and study of supersymmetric extensions of the 
SM with a stau NLSP. The rate of stau pairs has already been estimated in Refs. [7| 
and in]. In this paper, we will focus on the limited detector response for stau events using 
the detailed calculation of the stau energy loss in Ref. 0. We will carefully calculate 
the NLSP event rates for a cubic-kilometer neutrino observatory such as IceCube [TUj . 
taking into account the dependence of the detection efficiencies on stau energy and 

I The same is true in theories with an axino LSP, whose interactions are strongly suppressed by the 
large Peccei-Quinn scale, see e.g. p. The NLSP can also be very long-lived if its mass is very close 
to that of a neutralino LSP |2 • We do not study these alternatives in detail, but expect virtually the 
same results as in the case considered. 
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Figure 1. Chargino and neutralino exchange diagrams 
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spatial resolution. For the superpartner mass spectrum, we will use, for illustration, 
the benchmark point corresponding to SPS 7 from jTT] and a toy model where all 
superpartner masses are just above the experimental limits. Similarly, for the yet 
unknown high energy cosmic neutrino flux we will adopt as a benchmark the Waxman- 
Bahcall flux jl2| , which assumes a simultaneous production of cosmic neutrinos, protons, 
and gamma rays in astrophysical accelerator sources of the observed cosmic rays with 
comparable luminosities. 

The outline of the paper is as follows. In section |21 we calculate the differential 
production cross section of SUSY particles. We discuss their hfetime and energy loss 
in the Earth in section El Finally, we compute the resulting flux of staus and the event 
rates at neutrino observatories in section |3] 



2. NLSP Production 

Inelastic scattering on nucleons is the dominant interaction process of the high energy 
cosmic neutrinos in the atmosphere or the Earth. In this section, we calculate the 
differential cross sections for neutrino-quark scattering into sleptons and squarks. They 
determine the energy spectrum of the stau NLSPs. The resulting SUSY contribution to 
the neutrino-nucleon cross section is compared with the corresponding SM contribution. 
The parton-level SUSY contributions to charged and neutral current interactions 
between neutrinos and nucleons are the chargino and neutralino exchange diagrams 
shown in Fig. ^ as well as analogous interactions with the heavier quarks. The reactions 
produce sleptons and squarks, which promptly decay into the lighter stau, usually 
composed predominantly of the superpartner of the right-handed tau. The parton-level 
cross sections for these diagrams§ are given by 
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Tlie cross sections for tlie analogous processes withi anti-neutrinos are the same. 
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Figure 2. The SUSY mass spectrum of the benchmark point corresponding to 
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Summation over repeated indices is implied. We use the conventions and notations of 
|13| IT^ . The masses of the four neutrahnos and two charginos are denoted by m^o 
and m^-, respectively. The mixing matrices are Z^r and Z±. The neutralino couplings 
y = Y Z^sw + T^Zj^cw depend on the hypercharge Y = Q — T^ and the weak isospin 
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In the following we will focus on two SUSY mass spectra. One is given by the benchmark 
point corresponding to SPS 7 [TTj for a gauge-mediated SUSY breaking (GMSB) scheme 
with a messenger mass of 80 TeV. The corresponding mass spectrum calculated by 
SoftSusy 2.0.4 \TE\ is shown in Fig. |21 The other spectrum (denoted by "min m" 
in the following) consists of light charginos, neutralinos and sleptons at 100 GeV and 
squarks at 300 GeV. It is not motivated by any particular SUSY breaking scenario, but 
oriented at current experimental limits, in order to give an impression of what can be 
obtained in a very optimistic scenario. 
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Figure 3. The charged (red dotted) and neutral (blue dashed) current cross section 
for neutrino-nucleon scattering compared to the chargino and neutralino exchange for 
the SPS 7 benchmark point and a scenario with light squark masses of 300 GeV. 

The resulting neutrino-nucleon cross sections are shown in Fig. OJ compared to the 
SM contribution from W and Z exchange. Here and in the following computations we 
have used the CTEQ6D parton distribution functions ^H]- The contribution to the 
cross sections from charged currents is about twice as large as the one from neutral 
currents. 

In Eqs. ()la j) -()ld j) . we have neglected family mixings and contributions proportional 
to Yukawa couplings. We have taken into account the exchange of the heavier neutralinos 
and charginos, since generically their contribution to the cross sections is large and can 
even dominate. At the SPS 7 benchmark point we consider, assuming equal neutralino 
masses turns out to be a rather good approximation that helps simplify analytical 
calculations. Nevertheless, we use the exact expressions in the numerical calculations 
discussed in the following. 



3. NLSP Propagation 

The squarks and sleptons produced in the interactions discussed in the previous section 
promptly decay into stau NLSPs. Now, we discuss the survival probability of the latter 
after propagation towards the Cherenkov light detector. The corresponding effective 
detection volume is limited by energy loss and decay. 
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3.1. Energy Loss 

In the case of the muon, the mean energy loss per column depth 2;, measured in gcm~^, 
is given as 

^)^«, + /?M^M- (2) 

Here, a^ is determined by ionization effects and (3^ accounts for bremsstrahlung, pair- 
production, and photohadronic processes. In general, these coefficients are weakly 
energy dependent. For our purpose, it is sufficient to approximate their values by 
constants, a^ ~ 2 ■ 10^^ GeVcm^g"^ and /?^ ~ 4 ■ fO~^ cm^g^^. For muons above a 
critical energy E^^ = a^/P^ ~ 500 GeV, radiative energy losses dominate, which are 
proportional to the energy itself. Hence, Eq. (J2)) can be used to determine the muon 
energy if i?^ 3> E"^. We will take this critical energy as an effective cut-off in the 
following since the detection efficiency becomes very small below ^Tj. The range of a 
muon is then given as 

'^'^ P, \a, + P,E';i 

In the case of a stau, the radiative term in Eq. (J2I) is suppressed, since [3r satisfies 
PTfTTir ~ P^rn^. At large energies this gives an increase of the stau range compared to 
that of the muon by the mass ratio rrir/rn^. In the following we will use the results 
of [H] for the range of the stau.|| 

3.2. Decay 

The lifetime of the stau NLSP can be very long, since its decay into the gravitino LSP 
can proceed only gravitationally. The corresponding decay length L for relativistic staus, 
in units of the Earth's diameter 2/2®, is (e.g. [TH] ) 



R{E,,E-)^-\n{ / 7^; 1. (3) 



.2Re,) ^ VIOO GeV/ UoO keV/ V500 GeV/ 

The minimal stau energy we consider is about 500 GeV, corresponding to the critical 

energy of the muons. ^ Hence, staus with a mass not much larger than 100 GeV will 

always reach the detector, if the gravitino is heavier than 400 keV."*" In this case, the 

gravitino is a viable candidate for the dark matter [201 ■ Constraints from big bang 

nucleosynthesis and the cosmic microwave background yield an upper limit on the 

II If the stau NLSP has a large admixture from tl, the superpartner of the left-handed tau, charged 
and neutral current effects, which have not been considered in [H], might be important |18j . Assuming 
a small admixture of < 10%, as it is the case at the SPS 7 benchmark point, we expect that these 
effects can safely be neglected. 
% Note, however, that the critical energy of the stau is given as E" « {■mr/m^)E", according to 

+ Actually, we expect that our results remain unchanged even if the gravitino is lighter by an order of 
magnitude, since the dominant contribution to the event rate is due to staus with energies considerably 
larger than 500 GeV. 
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gravitino mass between 10 and 100 GeV [21]. Masses in the lower part of the allowed 
region are typical for models with gauge-mediated SUSY breaking, while masses of some 
tens of GeV can occur in gravity and gaugino mediation (for a review, cf. e.g. |22j). 

4. Detection Rate at Neutrino Telescopes 

There are several Cherenkov light based neutrino telescopes currently operating 
or under construction in the Antarctic (AMANDA P^l |21] and IceCube jTO]); in 
Lake Baikal [23 121, and in the Mediterranean (ANTARES [23 , NEMO ^, and 
NESTOR [2n]). The total flux of staus through a detector is proportional to the initial 
flux of high energy cosmic neutrinos F, attenuated inside the Earth according to the 
total inelastic neutrino-nucleon cross section (JuN- A slepton (squark) produced with an 
energy Ej (Eg) and a cross section cr^jv^'^ will promptly decay into a stau which looses 
energy by radiative processes according to Eq. ^ and reaches the detector with energy 
Er- The flux of staus per area A, time t, steradian Q, and energy Er is then given as 

Ztot -C/niax 

E ■ i=UQ 

Here, rrip is the proton mass, and the column depth z of the Earth is related to /, the 
distance to the detector in the line of sight, as —dz = p{l, 9)dL For the calculation 
of the total column depth of the Earth we use the density profile given in [SO] and a 
detector center at 1.9 km depth, as it is the case for IceCube. The measure \dEi/dEr\ 
takes into account the energy loss of a stau as well as the mean energy fraction it carries 
in the decay of the initial slepton or squark. For the SPS 7 benchmark point we have 
estimated {E^) ^ O.SEg and (i?f) ~ 0.5Ej, respectively. For the optimistic scenario 
with light sparticles this changes to (Er) ~ 0.5Eq and (Er) ~ Ej, respectively. 

For the flux of high energy cosmic neutrinos, we adopt the Waxman-Bahcall (WB) 
flux, ElF{E^) ^ 2 ■ 10"^ cm-2 s"^ sr"^ GeV per flavor |I2|. This flux is based on 
the assumption that the observed ultra high energy cosmic rays are protons from 
extragalactic astrophysical acceleration sites. In these violent environments protons, 
neutrinos, and gamma rays are produced with comparable luminosities. As a maximal 
neutrino energy we take lO^'' GeV. We have checked that an increase of this cut-off does 
not affect our results. 

The analysis of particle tracks observed in neutrino telescopes is calibrated for 
muons. Their energy Efj_ can be reconstructed by the energy loss ^E per column depth 
Az and Eq. 0. This is appropriate for muons with energies above the critical energy 
E'j^ = a^/Pfj_ ~ 500 GeV. A high-energy stau with energy E^f will deposit a much 
smaller energy fraction in the detector compared to a muon with the same energy, due 
to the reduced radiative energy loss. This has two effects. Firstly, the critical energy of 
a stau Ef = a^/ I3r ~ {m^/m^)E'^ is much higher compared to the muons. Secondly, 
due to the calibration of the detector, high-energy staus will be detected as muons with 
reduced energy E^ ^ {m^/m:f.)E:f. Note that this correctly identifies the boost factor of 
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Table 1. The expected number of muon and stau events for a one year observation 
at IceCube, assuming the Waxman-Bahcall neutrino flux per flavor, ElF{Ey) k, 
2 • 10~^ cm~2 s~i sr~^ GeV J2]- The effective areas of IceCube have been calculated 
by Monte-Carlo simulations and are given in jl7| . 



the particle 7 = Er/mr = E^/m^. Also the Cherenkov angle is consistently identified 
by cosOc oc (1 — 1/7^)"^/^. For that reason, we will present our results in Fig. |3 in 
terms of the detected energy E^. Of course, the total number of events does not change 
by this recalibration. Note that the critical energy for the stau has the same value for 
muons and staus in terms of the detected energy. For completeness, we show our results 
also in terms of the neutrino energy. 

4.1. One-Muon vs. One-Stau Events 

The upper panels of Fig. 0] show the differential fluxes of muons and staus through the 
IceCube detector. The calculation of the total rate of events requires the knowledge 
of the detection efficiency of the corresponding particle depending on its energy and 
direction. This is usually combined with the cross sectional area of the detector as 
an effective area A^q. The IceCube Collaboration has determined the effective area 
Agg(cos^, i?;j) for upward-going muons by Monte-Carlo simulations in ^^I- The angle- 
averaged values are shown in Tab. ^ together with the expected rates of muon events 
for the WB flux. The expected rates are in good agreement with those quoted in [T7j . 
if one takes into account the different normalizations of the initial neutrino fluxes. 

Despite the fact that the staus are produced in pairs, they may also produce one- 
stau events, for example if one particle of the pair misses the detector. These are 
indistinguishable from one-muon events, at least for high energies. A very conservative 
upper bound on the rate of one-stau events is given by the total rate of staus computed 
from Eq. ^ and a detection efficiency Al^i^Er) ~ A'^^^Ea). This upper bound is also 
shown in Tab.^and Fig.|3]for the two different SUSY mass spectra. It shows that in the 
fiducial energy region one-stau events are subdominant compared to the flux of muons, 
even for the most optimistic case of very light squarks with masses around 300 GeV. This 
seems to be generic, independent of the initial high-energy cosmic neutrino flux, as we 
have checked. Hence, the total rate of one-particle events can be used for reconstructing 
the neutrino flux without taking into account the contributions from NLSPs. This 
decouples the analysis of the extra-galactic neutrino flux from the search for parallel- 
stau events, to which we turn now. 
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Figure 4. Expected fluxes and event rates for the example of IceCube and the 
Waxman-Bahcall neutrino flux ^21- The results are shown in terms of the detected 
energy Ed (left panels) and the neutrino energy Ei, (right panels). We have used a 
SUSY mass spectrum with light squarks of 300 GeV (blue upper solid) as well as the 
one from the benchmark point corresponding to SPS 7 (green lower solid). 
Upper panels: Fluxes of upward-going muons (red dotted) and staus (sohd). For 
optimistic SUSY scenarios the spectrum of staus may dominate over the flux of muons 
at large neutrino energies (right panel). However, the detected spectrum of staus will 
be subdominant (left panel) due to the detector calibration explained in the text. 
Lower panels: Rates of parallel-stau events for an energy cut-off at Ef^^ = a^j (3^ 
(thick solid) and at E" = af/Pr = {mr/m^)E'r^ (thin solid), respectively. 



4-2. Rate of Parallel- Stau Events 

The detection efficiency of stau pairs, i.e. coincident parallel tracks, depends on the 
energies and directions of the staus, as in the case of single events, and also on their 
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separation. As a detailed Monte-Carlo simulation is beyond the scope of this study, we 
will use a simple approximation for the resulting angle-averaged effective area in the 
following. Firstly, we introduce a low energy cut-off for the detectable stau energy. As 
examples, we use the critical energies of the muon and of the stau. Secondly, we assume 
that the observatory can identify an event of coincident parallel stau tracks, if they are 
separated by more than 50 m jHI] and less than 1 km. Here, we approximate the opening 
angle by the angle between sleptons and squarks in the laboratory frame. If the pair 
event survives these cuts, we assume an effective area of 1 km^ for its detection. 

The lower panels of Fig. |3] show the expected rates of these pair events for the two 
SUSY scenarios. Since the staus of a pair have different energies, we show the average of 
the corresponding spectra. Note that a higher energy threshold in the detector does also 
affect the range of the second stau. This has an effect on the pair rates even if one stau 
is above the threshold as can be seen in the plots. The main reason for the suppression 
at large energies E^ is the requirement that the opening angle of the staus be compatible 
with the restriction on the separation. On the other hand, the distance of the interaction 
point from the detector has to be smaller than the corresponding range of the second 
stau. This mainly causes the reduction at small energies E^. For the Waxman-Bahcall 
neutrino flux and one year of observation at IceCube, we predict about 0.09 pair events 
for the benchmark point corresponding to SPS 7 and about 5 events for our toy model 
with light sparticle masses. This calculation assumes an effective low-energy cut-off at 
E'^/^ . We also show the effect of increasing the cut-off to Ef = {mr/m^)E'^^ . This might 
be necessary due to the reduced radiative energy loss of the stau. Our expectations are 
then reduced to 0.009 and 0.2 events, respectively. Currently, the IceCube Collaboration 
quotes an upper limit at 90% confidence level of E^F{E^) < 3 ■ 10~^ cm~^ s~^ sr~^ GeV 
per flavor for the diffuse flux of cosmic neutrinos [22] ■ This is approximately one order 
of magnitude above the Waxman-Bahcall flux taken for our calculation. If the cosmic 
neutrino flux does not lie much below this limit, the predicted event rates will increase 
by a factor 10, giving up to 50 pair events for the optimistic SUSY mass spectrum. 

The expected background of parallel muon pair events from random coincidences is 
bounded from above by the number of muons arriving within a (generous) time-window 
of 1 fis. Even this is several orders of magnitude below the stau pair event rate, 

, , < -mm 1 : 1 us ■ — -^ ~ C 10 . 6 

However, double muon events from processes like v^ + N — >■ ^~ + tt^ + X — >■ 
IJL~ + n'^ + 1/^ + X are expected to be more likely. We expect that these contributions 
are still subdominant compared to stau pairs due to meson-nucleon interactions in the 
Earth. This should be studied in more detail in the future. 

5. Summary and Discussion 

High energy cosmic neutrinos collide with matter in the Earth at center of mass energies 
beyond the capability of any earth-bound experiment. The attempt to measure the 
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cosmic ray fluxes in various observatories is therefore tiglitly connected to extrapolations 
of the SM interactions of these particles to very high energies. Besides, one can look for 
deviations from the SM predictions for the strength of the interactions or even for new 
particles. 

In this paper, we have examined the question whether neutrino telescopes could 
play a significant role in the discovery and study of supersymmetric extensions of 
the SM with a gravitino LSP and a long-lived stau NLSP. In these models, neutrino- 
nucleon interactions can produce pairs of staus which show up as parallel tracks in the 
detector 0. 

We have argued that events with a single stau in the detector are virtually 
indistinguishable from muon events at high energies. However, the reduced energy 
loss of staus in matter, which increases the effective volume, reduces also the detection 
efficiency in the telescope. As a result, single stau events play only a subdominant 
role compared to muon events. Consequently, the reconstruction of the initial high- 
energy neutrino flux from the total rate of muon-like events can proceed assuming SM 
interactions alone. 

If the superpartner masses are close to their experimental limits and if the cosmic 
neutrino flux is also close to its current experimental bound, we expect up to 50 
pair events per year in a cubic kilometer detector such as IceCube, with negligible 
background. Less favorable event rates are obtained for a less optimistic cut-off in the 
stau energy or for SUSY mass spectra of commonly used scenarios for SUSY breaking, 
mainly because the squarks are significantly heavier. 

If the gravitino mass is larger than 400 keV, as we have assumed, the stau decay 
length is larger than the Earth's diameter. Thus, the spectrum is independent of the 
gravitino mass. For lighter gravitinos, the rate of stau pairs will drop starting at the low- 
energy end of the spectrum. Therefore, it might be possible to obtain information about 
the gravitino mass, if the spectrum is measured very accurately and if the superpartner 
masses are known. However, with the small number of pair events we find, this appears 
challenging. 
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